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Abstract. The physical conditions in typical starburst galax- 
ies are investigated through critical infrared (IR) line ratios, 
as previously suggested by Lutz et al. (1998). The calcula- 



tions by a composite model which consistently accounts for the 
coupled effect of shock and photoionization by hot stars defi- 
nitely fit the observed line ratios of single objects and explain 
the observed relation between [O lv]/([Ne n]+0.44[Ne III]) and 
[Ne lll]/[Ne II]. The shock velocity and the gas density are the 
critical parameters. 

Most of the shocks are produced in low density-velocity 
(no ~ 100 cm~ 3 and V s ~ 50 — 100 kms" 1 ) clouds which 
represent the bulk of the ionized gas in starburst galaxies. How- 
ever, though they are by many orders less numerous, high- 
velocity (~ 400 - 600 kms" 1 ) shocks in dense (~ 500 - 800 
cm -3 ) clouds are necessary to reproduce the critical IR line 
ratios observed in the low-excitation Starburst Nucleus Galax- 
ies (SBNGs: M82, M83, NGC 253, NGC 3256, NGC 3690, 
and NGC 4945). These model predictions are in good agree- 
ment with the powerful starburst-driven superwinds and highly 
pressured ISM observed in SBNGs. On the contrary, the high- 
excitation Hll galaxies (II Zw 40 and NGC 5253) do not show 
any clear signature of large scale outflows of gas. This differ- 
ence between Hll galaxies and SBNGs can be interpreted in 
terms of temporal evolution of their starbursts. 

Key words: galaxies: individual - galaxies: ISM - galaxies: 
starburst - infrared: galaxies - line: formation - shock waves 



1. Introduction 

Detection of faint emission in the high excitation [Oiv] 
25.9/it n infrared (IR) line has been reported by Lutz et al. 
( 1998 ) in a number of well-known starburst galaxies from ob- 
servations with the Short Wavelength Spectrometer on board 
ISO. 

IR line ratios such as [Oiv] 25.9/itn / [Nell] 12.8/itn and 
[O IV] 25.9/im / [Ne III] 15.6/im were used previously (Lutz et 
al., 1996; Genzel et al., 1998) to establish the dominant source 



or [Ne III], Lutz et al. ( 1998 ) use [O lv]/([Ne II] + 0.44[Ne III]) 



as a critical line ratio to reveal different excitation conditions 
in starburst galaxies. Results from pure photoionization mod- 
els or from pure shock models are used by these authors in or- 
der to explain the observed line ratios. However, these models 
were not successful in reproducing the IR line ratios of sev- 
eral starburst galaxies. In fact, models with a power-law radia- 
tion producing enough [O IV] line give too high values for the 
[Ne m]/[Ne II] line ratio. Those with two-temperature (39000 
K + 80000 K) black body radiation fail in explaining the low- 
excitation starbursts. Finally, pure shock models could provide 
[O IV] line intensities similar to observations but not enough 
[Ne II] and [Ne III] which should come from the Hll region. 

These results indicate that a composite model, which ac- 
counts for the coupled effects of shocks and photoionization, 
must be used. Indeed, in a disturbed region such as that con- 
taining starbursts, it is expected that shocked gas is embedded 
in stellar radiation. This kind of composite model has success- 
fully explained the AGN spectral features (Viegas & Contini 
1997 and references therein) as well as those of single galaxies 
where the AGN and starburst activities are of comparable im- 
portance (e.g. the Circinus galaxy, Contini et al. 1998). These 
models were provided by the SUMA code (see Viegas & Con- 
tini 1994 and references therein). 

In the present study we suggest that shocks result from 
supernova explosions in the nuclear starburst region of the 
galaxy, which provide enough kinetic energy to create high ve- 
locity superwinds. These in turn will accelerate gas clouds on 
large scales. The physical conditions in the ionized interstellar 
clouds are determined by shocks coupled to photoionization by 
massive hot stars (O and Wolf-Rayet). 

This paper, rather than being an answer to Lutz et al. 



(starburst or AGN) of ionizing radiations in ultraluminous in- 
frared galaxies. Actually, rather than normalizing just to [Ne II] 



( Jl998| ), is aimed at carrying on their investigation. The distinc- 
tion between AG N and starburst act i vity in galaxies bein g al- 
ready established ( Lutz et al., 1996 ; |Genzel et al., 1998 ), our 
aim is to find out which models (provided by the SUMA code), 
better fit the observed trend between IR line ratios in starburst 
galaxies. We will not discuss in detail the single galaxy struc- 
tures, but will use the fit of calculated to observed line ratios 
to define the critical parameters (and their ranges of values) 
which prevail in the emitting gas. In other words, we will treat 
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Table 1. Parameters for single-cloud models 



Galaxy 


V s 


no 


U 


D 




(kms -1 ) 


(cm -3 ) 




(pc) 


II Zw 40 


50 


100 


0.02 


0.10 


NGC 5253 


100 


100 


0.01 


0.01 


NGC 3690 B/C 


300 


100 


0.01 


0.03 


NGC 3256 


400 


500 


0.01 


1.70 


M 82 


400 


500 


0.01 


2.70 


M83 


500 


500 


0.01 


0.70 


NGC 253 


600 


800 


0.01 


0.23 


NGC 4945 


600 


800 


0.01 


0.20 



systems of gaseous clouds moving radially outward from star- 
forming regions ionized by hot stars and shocks, avoiding a 
deep investigation on the internal structure of single objects. 

2. The results by a composite model 



As commented above, Lutz et al. (1998) concluded that shocks 
are the most plausible ionizing mechanism for explaining the 
[O IV] line, while the [Ne II] and [Ne III] come from the pho- 
toionized region. However, shock results should be obtained 
by consistent models concomitantly with photoionization from 
hot stars. The SUM A code consistently accounts for both the 
shock effect accompanying cloud expansion and photoioniza- 
tion by hot stars, in a plane-parallel symmetry. The geometrical 
thickness of the clouds range from 0.01 to 2.70 pc (see Table §, 
which is small compared with the radius of the starburst regions 
in galaxies. In our model, we assume that the gas clouds move 
radially outward from the starburst region, i.e. the hot star pop- 
ulation is facing the inner edges of the clouds while a shock 
develops on the opposite side. 

For all models an average stellar temperature, T*= 5 x 
10 4 K, for the hot star population, and a preshock magnetic 
field, Bo = 10 -4 gauss, are adopted, as well as cosmic abun- 
dances. Lower abundances should lead to a lower cooling rate 
of the gas downstream which could be reajusted by adopt- 
ing higher preshock densities. The other input parameters for 
SUMA: the shock velocity, V s , the preshock density, no, the 
ionization parameter^} U, and the geometrical thickness of the 
clouds, D, are chosen to vary within the ranges consistent with 
observations. 

2.1. Single-cloud models 

First we perform calculations by single-cloud models, which 
are compared with the observations of each galaxy in Fig. [j] 
The input parameters in the models are given in Table]]]. Filled 
squares represent the data from the Lutz et al. (1998) sample 
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Fig.l. [Oiv] 25.9/mi / ([Nell] 12.8/im + 0.44 [Nem] 
15.5/xm) versus the [Nelll] 15.5/mi / [Nell] 12.8/im line ra- 
tio. The observational data from Lutz et al. (1998) are indi- 
cated by filled squares. Open symbols correspond to results for 
single-cloud (squares, input parameters listed in Table [j]) and 
multiple-cloud (triangles, input parameters listed in Table |J) 
models. 



and open squares the results of model calculations. Upper lim- 
its are not included. As shown in Lutz et al. (1998, Fig. 2) they 
are located in the left region of Fig. |J The fit is very good for 
almost all galaxies. A perfect fit for individual data is senseless 
because many different conditions coexist in single objects. 
Notice indeed that NGC 4945 is a composite object with star- 
burst/Seyfert 2 nuclear activity. This could explain its higher 
intensity of [Oiv] 25.9/im relative to model predictions for 
starburst only. A possible hidden AGN has also been reported 
in the merger galaxy NGC 3256 ( Kotilainen et al., 1996[ ). 

The parameters listed in Table y indicate that the shock ve- 
locity and the preshock density are the critical parameters to 
reproduce the observed trend of the IR line ratios in Fig. [j] 
The velocity growth is generally followed by an increase of 
the preshock density. An ionization parameter of 0.01 was 
adopted for nearly all the models except for II Zw 40 for which 
U = 0.02. For this galaxy the effect of radiation dominates 
the shock effect (see Sect. Q). The geometrical thickness of the 
clouds, D, ranges between 0.01 pc for NGC 5253 and < 3 pc 
for M 82. 



The ionization parameter is the ratio of the photon density of the 
ionizing radiations reaching the nebula to the gas density in the inner 



Lutz et al. ( 1998 ) already suggested that shocks are the 
most plausible source of [Oiv] emission observed in star- 
burst galaxies. Our results show that shocks are also neces- 
sary to reproduce the low observed values of [Ne m]/[Ne II] 
(see Sect. [| for details). Moreover, Fig. [l] shows that the ob- 
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Table 2. Calculated line fluxes (in erg cm 2 s 1 ) with multi-cloud models 



Line 






M 82 








NGC 3256 










NGC 253 






(Mm) 


Al 


A2 


A3 


A4 


Av 


Bl 


Bl 


B3 


Av 


CI 




C2 




C3 


Av 


[Cll] 158 


0.05 


0.026 


0.007 


0.420 


2.56 


0.03 


0.012 


0.007 


1.6 


0.03 




0.07 


6 x 


io- 4 


7.06 




0.01 


16.5 


27.0 


56.5 




0.04 


75 


24.5 




0.02 




99.1 




0.85 




[Om] 88 


9.6 


0.002 


0.008 


0.001 


0.92 


6.2 


0.002 


0.008 


0.8 


2.10 


3 X 


10 -5 




0.0076 


0.78 




6.30 


4.2 


89.0 


0.39 




13.9 


23.8 


62.4 




11.3 




0.3 




88.5 




[Oil 63 


0.004 


0.067 


0.007 


0.348 


2.96 


0.002 


0.025 


0.007 


2.9 


0.78 




0.075 


5 X 


10~ 4 


7.59 




0.0008 


37.0 


23.0 


40.0 




0.002 


86.5 


13.5 




0.5 




98.8 




0.66 




[Om] 52 


89.6 


0.004 


0.008 


0.003 


1.46 


56.6 


0.004 


0.008 


2.0 


20 


6.7 X 


1Q — 5 




0.0078 


1.73 




36.7 


43 


57.8 


0.9 




54.6 


20.4 


26.5 




53.7 




0.4 




45.7 




[Si ii] 35 


25.6 


0.166 


0.027 


2.170 


13.20 


15.9 


0.073 


0.027 


9.4 


3.70 




0.44 




0.01 


45.2 




1.2 


20.8 


20.5 


57.6 




3.2 


77.4 


19.4 




0.4 




97.4 




2.2 




[S in] 33 


12.0 


0.055 


0.001 


0.570 


3.10 


7.5 


0.024 


0.001 


2.6 


3.30 




0.02 


2 X 


IO" 4 


2.17 




2.3 


29.3 


4.2 


64 




5.5 


91.6 


3.0 




7.1 




91.9 




0.92 




[Sin] 18 


96.0 


0.031 


6.6 X IO -4 


0.114 


1.54 


60.6 


0.013 6.6 x 


10~ 4 


2.5 


30 




0.005 


1.2 X 


10~ 5 


1.85 




1.9 


32.9 


4.2 


25.8 




46.2 


52.2 


0.02 




75 




24.5 




0.64 




[Oiv] 25.9 


7 


1.2 X 10~ 4 


6.6 X 10~ 5 


2.3 x 10 ~ 4 


0.05 


7 


1.2 x 10 -4 6.6 X 


10 -5 


0.14 


5.6 


8.8 X 


10 -5 


6.3 X 


10 -5 


0.3 




82 


4 


13 


1 




89 


8 


3 




94 




3 








[Nem] 15.5 


123 


0.014 


0.008 


0.01 


1.80 


80 


0.0075 


0.008 


2.77 


46 


8 X 


io- 4 




0.0063 


2.8 




40 


13 


45 






55 


27 


18 




75 




3 








[Neil] 12.8 


960 


0.025 


0.004 


0.15 


7.06 


606 


0.01 


0.004 


12.6 


490 




0.0055 




0.0016 


23.3 




82 


6 


5 


7 




90 


8 


2 




97 








0.7 




W (%) 


0.005 


13.8 


83.3 


2.895 




0.01 


60.4 


39.59 




0.007 




49.493 




50.5 




V s (kms- 1 ) 


400 


200 


100 


200 




400 


200 


100 




600 




100 




100 




»0<cm-3) 


500 


50 


40 


100 




500 


50 


40 




800 




200 




40 




U 


0.01 


0.01 


0.01 


0.01 




0.01 


0.01 


0.01 




0.01 




0.001 




0.01 




D(pc) 


2.67 


3.33 


3.33 


20 




1.67 


1.67 


3.33 




0.23 




10 




0.5 





For each line, we report as a second entry the percentage contribution of each model in the weighted sum flux. 



Table 3. Observed vs. calculated IR line fluxes 3 



Line 


NGC 253 




NGC 3256 


M 82 




(A™) 


Obs b 


Cal 


Obs b 


Cal 


Obs b 


Cal 


[Oiii]88 


0.10±0.02 


0.1 


0.4±0.3 


0.5 


0.6±0.2 


0.4 


[O i]63 


0.9±0.1 


1.0 


1.2±0.4 


1.7 


1.1 ±0.4 


1.2 


[Om]52 


0.20±0.04 


0.2 


<1.3 


1.2 


0.7±0.2 


0.5 


[Sin]35 


1.1±0.2 


6.2 


1.9±0.5 


5.6 


1.0±0.4 


5.4 


[S m]33 


<1.3 


0.3 


1.2±0.3 


1.6 


1.1 ±0.3 


1.2 


[S in] 18 




0.2 


0.9±0.2 


1.5 


0.7±0.3 


0.6 



a Fluxes are relative to [C II] 158^im. 

b Observed line fluxes from Carralet al. ( |l994| ) for NGC 253 and 
NGC 3256, Rigopoulou et al. ( |l996| ) for NGC 3256, and Lord et al. 
( |1996| ) for M82. 



acterized by high velocity and density clouds. However, pre- 
vious investigations on the physical conditions in the ISM of 
NGC 253 (|Carral et al. 1994] ), NG C 3256 (parral et al., 1994] ; 
Rigopoulou et al., 1996| ), andM 82 ( |Houck et al., 1984 Lordet 
al., 1996) found relatively low electronic densities on the basis 
of the [S III] 33/im/18/im and [Om] 88^m/52/im line ratios, 
and low velocities derived from the FWHM of the IR line pro- 
files. 



2.2. Multi-cloud models 

Many different conditions coexist in the ISM of single galax- 
ies. Therefore, we performed multi-cloud model calculations 
to reproduce all the observed IR line ratios available in the lit- 



82. In these models, we assume that both low and high density- 
velocity (d-v) clouds are present in single objects. 

In Table ||, the calculated IR line fluxes are given for single- 
cloud models and for the weighted sums of the multi-cloud 
models. Notice that the line fluxes are calculated at the emit- 
ting nebula while the observations are given at Earth. The input 
parameters of single-cloud models are also listed. 

Line fluxes are proportional to the square of the gas emit- 
ting density, and compression downstream increases with the 
shock velocity. Thus, the line fluxes computed from high d-v 
models (i.e. Al, Bl, and CI) are higher than those calculated 
from low d-v models (A2, A3, A4, B2, B3, C2, and C3) by 
some orders of magnitude, except for lines [O I] and [C II]. The 
gas in high d-v models is ionized to very high temperatures 
(Fig. ||b), and so low-level lines like [O I] and [C II] are weak. 
Recall that the first ionization potential of C is lower than that 
of H (and Ne). Consequently, [C II] lines, as well as [O I] lines, 
are emitted by relatively cold gas and can be weak in models 
which show strong [Ne II] lines. 

In order to fit the observations, a balance must be achieved 
between high and low d-v models in the multi-cloud spectra, 
and the final result is a weighted sum (A v in Table 2) with the 
weights chosen adequately. For each model, the weight W rep- 
resents the relative number of clouds (in %) characterized by 
that model. The weights of low d-v models must be by a factor 
of 10 4 — 10 5 higher than for high d-v models. In other words, 
the number of high d-v clouds must be very low compared with 
the number of low d-v clouds. In Table ||, for each emission 
line, the second entry corresponds to the percentage contribu- 
tion of the single-cloud model to the line flux, taking into ac- 
count the corresponding weight. The high d-v clouds contribute 
up to 90% to the [O IV] and [Ne II] emission lines, whereas the 
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d-v clouds for M 82 and NGC 3256. The others IR lines are 
mainly emitted by the low d-v clouds except the [O III] 52 line. 

The results of multi-cloud model calculations are compared 
to the observed IR line ratios in Table ^| and plotted in Fig. [l] 
(open triangles). The fit is good for all line ratios (especially for 
the density critical ones, see sect. |2.1[ ), except for the [Si II] line 
which is overpredicted by a factor of about 5. This might indi- 
cate that silicium is depleted and locked in dust grains. Adopt- 
ing a relative abundance Si/H lower by a factor of ~ 5 than cos- 
mic will not affect the predictions for the other line ratios listed 
in Table |. Most of the IR lines reported in Table | (e.g. [O I], 
[Cll], etc) are mainly produced by the low d-v clouds. Thus, 
in the weighted sum, the IR lines ratios [S III] 33/im/18/im and 
[Olll] 88/im/52/im are determined by low d-v models, while 
the [Olv]/([Nell]+0.44[Nelll]) and [Ne lll]/[Ne II] line ratios 
are mainly determined by the high d-v cloud component. 

Once fitted the line ratios, it is now possible to examine the 
FWHM of the line profiles, because the weights adopted in Ta- 
ble H for high d-v models are very low and could lead to hardly 
observable broad components in the line profiles. The calcu- 
lated line profiles for M 82 and NGC 253 are shown in Fig ||. 
The effect of galaxy rotation is included in the calculation. Ro- 
tation velocities are about 76 km s _1 and 189 km s _1 for M 82 
and NGC 253 respectively. Two significant lines are chosen for 
each galaxy: [O IV] and [Ne III] for M 82, and [O IV] and [O III] 
88 for NGC 253. The [O III] 88 line is chosen because it repre- 
sents the case of a low contribution of the high d-v model (CI) 
relatively to a strong contribution of the low d-v model (C3). 
Single-cloud model profiles are also plotted in the figures us- 
ing their percentage contribution reported in Table |[ The broad 
component from the high d-v clouds dominates the weighted 
averaged profile of [Oiv], while in the [Olll] 88 and [Nelll] 
profiles the broad component is faint and hardly observable, 
due to the low contribution from high d-v clouds. For M 82, 
the predicted FWHM of the [Oiv] profile is - 300 kms~\ 
in good agreement with the ob serve d value (FWHM ~ 360 
kms -1 ) reported by Lutz et al. (1998). 



3. The physical conditions across emitting clouds 

To better understand the role of the input parameters in deter- 
mining the physical conditions of the emitting gas, the distri- 
bution of some important physical quantities across the cloud 
is given in Figures |]a and ||b respectively for two significant 
models listed in Table [if NGC 5253, which represents the 
low density-velocity cloud component, and the high density- 
velocity component of NGC 253. Each figure shows three pan- 
els where the shock front is on the left and the photoionized 
edge is on the right. The top panel displays the electron tem- 
perature, T e , and density, n c . The fractional abundances of the 
ions +3 /0, Ne +2 /Ne, and Ne +1 /Ne are shown in the central 
panel. The ionization rates for Ne +2 due to the primary radia- 
tion (P), to the diffuse radiation (D), and to electronic collisions 
(C) are given in the bottom panel. 

The [O IV] line is always emitted from the shock dominated 



NGC 253 
[OIV] 25-9 




-500 
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500 


i 1 i i i 




i i | i i i | i 
M82 






[Nelll] 15.5 _ 









Fig. 2. Predicted line profiles for NGC 253 (top) and M 82 (bot- 
tom) using multi-cloud models. The final profile (solid line) is 
the weighted sum of different cloud components: high density- 
velocity clouds (dotted line) and low density-velocity clouds 
(long-dashed line: V s = 200 kms -1 ; short-dashed line: V s = 
100 kms" 1 ) 

from the photoionized zone. In particular, Fig. |^a shows that 
the critical temperatures are about 10 5 K. For NGC 5253 the 
shock effect becomes negligible at a distance of about 10 15 
cm from the shock front, and most of the physical conditions 
in the cloud are then dominated by the effect of stellar radia- 
tion which leads to [Ne m]/[Ne II] > 1 . On the other hand, due 
to strong compression in the NGC 253 gas clouds (Fig. ^]b), 
the temperature rapidly drops to T c < 10 4 K. This can be ex- 
plained by the high density (~ 10 5 cm~ 3 ) of the gas in the 
innermost region of the cloud which strengthens the cooling 
rate. The role of the primary radiation is to keep the gas at a 
temperature slightly below 10 4 K corresponding to a relatively 
high Ne +1 /Ne compared to Ne +0 /Ne and Ne +2 /Ne. Therefore, 
for this galaxy, [Ne m]/[Ne II] is < 1 . Owing to the high opti- 
cal thickness of the gas, the stellar radiation flux reaching the 
inner edge of the cloud rapidly decreases, the size of the Ne +2 
photoionized zone is reduced, and the diffuse radiation effect 
prevails accross the cloud. 

For NGC 253, temperatures as high as 7 x 10 6 K deduced 
from the ROSAT PSPC data (Heckman 1996 and references 
therein) are in agreement with 5.4 x 10 6 K obtained by the 
model in the postshock region (Fig. ||b). Precursor soft X-rays 
radiation from the hot gas could affect the ionization conditions 
in the upstream gas clouds. However, due to the adiabatic jump 
and compression, the densities n upstream are lower than in the 
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Fig. 3. The physical conditions across low-density (a, NGC 5253) and high-density (b, NGC 253) clouds are shown in three 
panels, each one with a vertical line indicating the center of the emitting clouds. The electron temperature (dotted line) and 
electron density (solid line) appear in the top panel. The fractional abundances of +3 (long-dashed line), Ne +2 (solid line), and 
Nc +1 (short-dashed line) are displayed in the central panel. The ionization rates for Ne +2 (in logarithmic units) due to electron 
collisions C (long-dashed line), to the stellar radiation P (solid line) and to the diffuse radiation D (short-dashed line) are shown 
in the bottom panel. The shock front is on the left side and the edge of the cloud photoionized by the stellar radiation is on the 
right of each panel. The horizontal axis scale is logarithmic and symmetric in order to show an equal view of both sides of the 
emitting cloud 



intensities depend on n 2 , the contribution of the line emission 
from upstream gas to the final spectrum will be small. 

4. Comparison with the observed starburst properties 

One might interpret the observed trend in Fig. [l] as a continuous 
sequence between two extreme cases: starburst regions where 
shocks in high-velocity clouds are necessary to reproduce the 
observed critical IR line ratios (bottom left part of the diagram) 
and the "radiation-dominated" ones (upper right part of the di- 
agram), where the shock contribution is relatively less impor- 
tant, because formed only in low-velocity clouds. Let us now 
check wether the critical parameters derived from the models 
along this sequence are in agreement with the observed prop- 
erties of the starburst galaxies. 

Galaxies in the first group (the last six galaxies in Table Q) 
are massive and chemically evolved spiral galaxies which ex- 
hibit a single energetic and compact starburst in their nucleus, 
and are commonly named Starburst Nucleus Galaxies (SBNGs; 
Balzano 1983; Coziol et al. 1994). Evidence for powerful large- 
scale starburst-driven superwinds in the nuclear region of SB- 
NGs (like M 82, NGC 253, NGC 4945 and NGC 3690 in the 



Table 4. Galaxian properties 



present sample) has been reported by Heckman et al. Q1990J) 
Moreover, the strong X-ray emission (Lx ~ 10 40 
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Galaxy 


Morph. a 




d b 
(Mpc) 


(ergs 1 ) 


II Zw 40 


BCDG 


-14.5 


9.2 




NGC 5253 


Irr 


-16.2 


4.1 


38.6 


NGC 3690 B/C 


Merger 


-21.1 


42.4 


41.0 


NGC 3256 


Merger 


-22.2 


33.2 


41.7 


M82 


Sd 


-19.0 


4.3 


39.7 


M83 


SABc 


-21.2 


5.1 


40.2 


NGC 253 


SABc 


-20.9 


3.3 


39.9 


NGC 4945 


SBcd 


-20.1 


5.1 


39.9 



see 



a From the LEDA database 

b The distance (d) is derived from the Galactic Standard of Rest 
(GSR) velocity, using a Hubble constant Ho — 75 km s" 1 Mpc -1 , 
except for the distance of NGC 5253 taken from Saha et al. (1995). 
The X-ray luminosity is derived from the flux reported in the ROSAT 
catalog (Brinkmann et al. 1994) or in the EINSTEIN catalog 
(Fabbiano et al. 1992). 



of hot superbubbles due to the combined action of stellar winds 
from massive stars and supernova explosions in the central star- 
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km s 1 have been estimated for the ionized gas in M 82 (Shop- 
bell & Bland-Hawtho rn, 1998), NGC 4945 (Moorwood et al., 
1996), and NGC 253 ( fceckman, 1996| ). 

In the central region of SBNGs, part of the ISM is over- 
pressured relative to the ISM in the Milky Way (Heckman et 
al., 1990; Carral et al., 1994), due to compression by super- 



winds. The nuclear starbursts usually contain many compact 
radio sources associated with bright and small supernova rem- 
nants (Muxlow et al. 1994; Huang et al. 1994 for M 82; Ul- 
vestad & Antonucci 1997 for NGC 253; Zhao et al. 1997 for 
NGC 3690). Such a population of dense and compact objects 
is generally attributed to the high pressure in the central region 
of SBNGs which allows the formation of dense ionized clouds 
(see e.g. Carilli 1996). 

The galaxies in the second group (the first two galaxies in 
Table are mostly metal-poor dwarf irregular galaxies with 
multiple star-forming r egions distrib uted randomly, and are 
known as Hll galaxies ( French, 1980 ). The archetypes of this 
class of galaxies are II Zw 40 and NGC 5253. Their starbursts 
are very young (< 4 Myr), as proved by the numerous Wolf- 



Rayet stars still present in II Zw 40 (Kunth & Sargent, 1981 
Vacca & Conti, 19921) and NGC 5253 dSchaerer et al.. 



1997 



Schaerer et al., 1999), and only contain a few supernova rem- 



nants (Vanzi et al. 1996 for II Zw 40; Beck et al. 1996 for NGC 
5253). Moreover, there is no evidence for powerful starburst- 
driven superwinds in these galaxies (Martin 1998 for II Zw 40; 
Marlowe et al. 1995 for NGC 5253) which are faint or even 
undetected in X-rays (see Table Q). 

Thus, the distinction between Hll galaxies and SBNGs can 
be understood in terms of temporal evolution of the starbursts. 
In the former galaxies (II Zw 40 and NGC 5253), the star- 
forming event is still young enough to contain a high propor- 
tion of massive O and Wolf-Rayet stars. These populations of 
hot (T c ff ~ 60000 K) stars produce large amounts of ioniz- 
ing radiation, but did not have the time to evolve off the main 
sequence to produce supernova explosions. In SBNGs, like M 
82 or NGC 253, the starburst is much older (> 20 Myr, e.g. 
Engelbracht et al. 1998), and the most massive stars had the 
time to evolve through supern ovae. As shown by the models of 
Leitherer & Heckman ( 1995 ), the ratio between the mechani- 
cal energy and the ionizing Lyman continuum, both injected in 
the ISM by the massive star population, rapidly increases by 
a factor of ~ 100 between ~ 5 and 20 Myr. This explain the 
difference between the "radiation-dominated" Hll galaxies and 
the "shock-dominated" SBNGs. The presence of high density- 
velocity clouds in SBNGs, which host a relatively "old" star- 
burst, require large quantities of mechanical energy to power 
large-scale superwinds. These energetic events disrupt molec- 
ular clouds through shock waves and drive them out of the nu- 
cleus at high velocities. 

5. Conclusions 

The observed relation between [O lv]/([Ne ll]+ 0.44[Nelll]) 
and [Ne lll]/[Ne II], and single object IR line ratios are well 



and photoionization effects occurring in starburst regions. The 
critical parameters which reproduce the observed trend of the 
IR line ratios are the preshock density and the shock veloc- 
ity which are consistently correlated. Most of the shocks are 
produced in low density- velocity (no ~ 100 cm -3 and V s 
~ 50 — 100 kms' 1 ) clouds which represent the bulk of the 
ionized gas in starburst galaxies. However, even if they are 
by many orders less numerous, high-velocity (~ 400 — 600 
kms -1 ) shocks in dense (~ 500 — 800 cm~ 3 ) clouds are nec- 
essary to reproduce the critical IR line ratios observed in the 
low-excitation SBNGs. 

These model predictions are in good agreement with the 
observed ISM properties in starburst galaxies. SBNGs are the 
only ones to exhibit powerful starburst-driven superwinds and 
highly pressured ISM. On the contrary, the high-excitation Hll 
galaxies do not show any clear signature of large scale outflows 
of gas. This difference between Hll galaxies and SBNGs can be 
interpreted in terms of temporal evolution of their starbursts. In 
the former, the starburst is too young to release large amounts 
of mechanical energy through supernova explosions whereas in 
the latter, characterized by older starbursts, the kinetic energy 
compete with the ionizing radiation in the heating of the ISM. 
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